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Abstract 

The behavior of lipid membranes in the presence of an external electric field is studied and used 
to examine the influence of such fields on membrane parameters such as roughness and show that 
for a micro sized membrane, roughness grows as the field increases. The dependence of bending 
rigidity on the electric field is also studied and an estimation of thickness of the accumulated 
charges around lipid membranes in a free-salt solution is presented. 
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1 Introduction 

Membranes can be modeled as statistical-mechanical systems, showing a variety of configurations 
as 2D fluid-like surfaces. In addition to thermal fluctuations, the fluctuation spectrum of a lipid 
membrane can be influenced by several other external factors such as osmotic pressure or pH dif- 
ferences. Another significant factor contributing to these fluctuations is the presence of an external 
electric field which can induce dramatic changes on the way a membrane behaves, rupture or shape 
transformation [1] being one example. This is due to the fact that when an electric field is applied 
to a lipid membrane it causes an extra lateral tension (referred to as electrotension) to appear which 
depends on the strength of the applied field as well as on a few constitutive parameters of the lipid 
membrane [2, 3, 1, 4]. 

Electrostatic interactions influence membrane elastic parameters in a fundamental way, in par- 
ticular the stability of flexible membranes [5] . The non- vanishing excess membrane charge causes an 
increase in membrane undulations due to Coulomb repulsion, while charge fluctuations or free ions 
in a solution with screening effects suppress this instability. It has been shown that [5] with respect 
to free- ion or ion solvents, the surface tension always decreases whereas bending rigidity may either 
decrease or increase. Such instabilities are due to changes in elastic parameters and may lead to a 
roughening of the lipid membrane surface when exposed to thermal and electrical fluctuations. 

Roughening effects have important consequences on mechanisms governing the motion of molecules 
or particles along the bio-membrane surface. Membrane roughness affects the interaction between 
the surface and outer particles in one hand and the motion of external inclusions [6, 7] on the other. 
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Also, roughness fluctuations affect the electronic properties of a membrane. These fluctuations are 
observed in graphene sheet and have an important role in its electronic properties [8]. In this regard, 
the relevant parameter attributed to a membrane is the roughness, defined as the root-mean-square 
value of vertical variations [6] and depends on several membrane parameters such as bending rigidity, 
surface tension, linear size, temperature and so on [9, 10]. For a detailed discussion sec [7] and the 
references therein. Molecular and monte carlo simulations have also been used to extract Nano- 
Scale lipid membrane roughness using various kinds of interaction parameters [11, 12]. Apart from 
roughness, other phenomenon such as electrop oration in the case of giant vesicles or a patch of 
membranes [13], lead to important biotechnology applications. 

Charge currents across a nearly flat and poorly conductive lipid membrane give rise to out-of- 
equilibrium membrane undulations. In this regard, two variant membranes under the influence of 
an electric field, namely, a fluid membrane attached to a rigid frame and a freely floating membrane 
can be considered. The theory behind the above mentioned cases may be inferred from a general 
situation where an infinite membrane is subjected to an external electric field. In such 
effective negative surface tension appears in the membrane, making it tense yet floppy-looking [14]. 
This theory is our preferred approach in computing and describing the roughening and other quantities 
relating to lipid membranes. The mean dynamical roughness of the membrane can be obtained by 
calculating the height-height correlation function [15]. 

In previous works [16, 17], we computed the stochastic trajectories of objects (inclusions) existing 
within and on the surface of a membrane via the application of Langevin dynamics. Also, by intro- 
ducing the height function as a stochastic Wiener process, a relation between random fluctuations of 
height and lateral diffusion of membranes was studied [17]. In this paper, we investigate the effects 
of electric fields on surface tension and bending rigidity and their influence on the roughness in lipid 
membranes. Our analysis is based on the extension of the results previously presented in [5] and [14] 
for the bending rigidity and surface tension. For simplicity, we employ an almost infinite flat lipid 
membrane subject to an electric field. One result is that the weaker the electric field, the smaller 
the roughness. If a slab made of charges with a certain thickness is considered around the undulated 
membrane, the dependence of the thickness on the electric field and other lipid membrane parameters 
can be estimated. 

The paper is organized as follows: in section two we write the modified Canham-Helfrich Hamil- 
tonian where the surface tension and bending rigidity include a term resulting from the application of 
the electric field. This Hamiltonian is then used to define the dynamical roughness of the membrane 
in section three and conclusions are drawn in the last section. 

2 Energetics of an almost plane membrane 

The free elastic energy of a symmetric, nearly flat, membrane is described by the Canham-Helfrich 
Hamiltonian [18, 19]. To study a nearly flat membrane, it is convenient to consider it parallel to 
the (xi,X2) plane, regarded as the reference plane. A single-valued height function h represents the 
position of a point on a fluctuating, nearly flat, sheet relative to the reference plane and in this, the 
so-called Monge representation [15] the Hamiltonian is written as 



where cto is the surface tension and hq is the bending rigidity of the membrane. In this form, 
the Hamiltonian is expressed solely in terms of the height function, h, and its derivatives. The 
Hamiltonian in (1) describes the energetics of the membrane from which one may obtain the membrane 
roughness. Since our main goal in this work is to study the behavior of lipid membranes in the presence 
of an external electric field, it would be advantageous to write the above Hamiltonian when such a 
field is present. Therefore, in what follows, we study a membrane under the influence of an electric 
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Figure 1: Net accumulation of conduction charges near a curved lipid membrane in an electric field, 
figure is based on ref [14]. T e i is the resulting normal stress component in the z-direction which 
enhances the local undulations h(r) of the membrane. It also induce a further surface tension in the 
membrane surface. 

field and investigate the resulting effects that such a field may have, using a modified form of the 
above Hamiltonian. 

2.1 The effective surface tension in an electric field 

Since lipid bilayers are impermeable to ions, in the presence of an external electric field, charges 
accumulate at the bilayer interface, as shown in Fig. 1. Thus, an additional transmembrane potential 
is created across the membrane. An almost flat, infinite and poorly conductive membrane under the 
influence of a fixed external electric field may be used as a simple model to investigate the effects of 
the electric field on lipid membranes, see Fig. 1. 

One may embark on obtaining a relation between the local undulation of lipid membranes in the 
large-wavelength limit by using the discontinuity of the Maxwell stress tensor across the membrane 
interface and solving the corresponding electro-kinetic problem by employing the usual boundary 
conditions on displacement vector and current density. This method was employed in [14], leading to 
a net decrease in energy at long wavelength limit relative to the membrane thickness d. This extra 
electric field contribution to the Canham-Helfrich Hamiltonian can be written as 

n el = ^J Sx{vhf . (2) 

Here a e i is an induced effective negative surface tension 

a el = -aE 2 , (3) 

where 
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a = d e ri . 

-9m . 

and E is the electric field strength. As for the other quantity e m refers to the permittivity of the 
membrane and g m , g\ are the conductivities of the lipid membrane and the surrounding medium 
respectively. 

The similarity of the equation (2) and the second term on the RHS of equation (1) is clear. Both 
represent the surface tension contributions to the Hamiltonian, one relating to the cytoskeleton- 
generated surface tension and the other only appears when an external electric field is applied. 
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Figure 2: Total surface tension versus the electric field for three different values of oq. 

For typical values of the parameters quoted in [14] , one gets q = 0.71 x 10~ 9 ^. Therefore, the 
total membrane surface tension at (which appear in the surface tension part of the Hamiltonian) as 
the sum of electricaly induced part a e \ and cytoskeleton-generated membrane surface tension ao may 
be written as 

ft = <?el + CO- (4) 

The electric part of surface tension a e i, being negative, enhances the undulations of the membrane 
while o"o, the intramembrane tension, tends to keep the membrane surface unaltered. Note that we 
use positive values for at to satisfy the nearly-flat assumption. Typical values for the intramembrane 
tension and bending rigidity of a fluctuating membrane at room temperature for different lipid ma- 
terials are 10" 9 ^ < a < W~ 6 ^i [17] and k ~ 10- 20 J (~ 2.5k B T) [20], respectively. These values 
correspond to the following range of values for the electric field, using equation (3) and (4), presented 
in Table I. 
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Table I: Upper bounds for the electric field corresponding to arbitrary intramembrane surface tension ao, used in the 
present work. As can be seen, the weaker the electric field, the lower the value of ao. A membrane with larger ao can 

withstand higher electric fields. 

The above values for the electric field have been taken with the view that the assumption of a nearly 
flat membrane is satisfied and is not violated by the application of too strong an electric field. They 
indicate that a membrane with weak surface tension is more affected by an electric field than a 
membrane with a stronger surface tension. Figure 2 shows variations of at with respect to the electric 
field for three values of a^. 

2.2 The effective bending rigidity in an electric field 

One can assume a similar equation for the total bending rigidity Kt as the sum of the electric contri- 
bution K e i and the cytoskeleton-generated rigidity kq, that is 

Kt = K e l + K . (5) 

In what follows, we will make an estimate of K e i using the theory of charged membranes and equation 
(4) and show that K e \ is two order of magnitude smaller than kq used below. 
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Let us assume that the membrane is charged due to the influence of an external electric field and 
the resulting accumulated net charge fluctuates on the surface. In such a scenario the contributions 
of a e i and k e i to the surface tension and bending rigidity which stem from charge fluctuations, are 
both negative. In the limit of large membrane undulations, the following relation for the negative 
surface tension may be written [5] 



2ttA , / 2vrA 
In 1 + 



(6) 



where A is a characteristic length similar to that of the Gouy-Chapman length especially in the 
presence of counterions around the membrane and a is the molecular cut-off. The solution of the 
Poisson-Boltzmann equation for a simple charged surface in the presence of counterions would enable 
one to define a characteristic length which decreases when the surface charge density goes to higher 
values. In the present study the membrane is charged via the external electric field and is not flat. 
However, we have assumed that the situation is similar to the case of a charged surface when A S> a. 
In this limit, we may safely ignore the logarithm term in equation (6) against 27rA/a. Equations (3) 
and (6) then implicitly lead to an expression for the characteristic length, A, mentioned earlier 

A^f^V. (7) 

As can be seen, an increasing electric field causes A to decrease as a result of the increasing charge 
density, leading to a highly charged membrane. Also, when a — > 0, one has A — > oo, showing that 
for a small molecular cut-off, the membrane becomes poorly charged and that in the presence of 
counterions, the charges do not accumulate near the membrane and often move freely in the solvent. 
On the other hand, if ^ decreases or n m increases, A increases, since there would be a concentration 
of charges on the boundaries. The direct dependence of A on temperature is clearly seen. The 
dependence of A on the electric field strength and membrane conductivity is shown in Fig. 3. 

In the same manner, if we assume that the membrane is charged by the application of the electric 
field so that net charges on the membrane can fluctuate, in the long wavelength limit, the following 
expression for the reduced bending rigidity holds [5] 

where L is a typical membrane length taken as 40 /im here. Figure 4 shows variation of K e \ with the 
electric field and g rn for g = 1^, a = 2nm, T = 300K and e m = 2eo- Note that the membranes 
become more flexible for higher values of both the field strength and membrane conductivity. This 
flexibility decreases rapidly when the membrane conductivity is low and the electric field strength is 
weak. It should be emphasized that the theory explained above may also be used when other sources 
are present. It has been used to study membranes with excess surface charge density [5] by dividing 
the Hamiltonian into three terms, that is, electrostatic interaction, entropic contribution and elastic 
energy. The effects of ionic salts on bending rigidity around cell membranes and in particular of 
micron-sized vesicles are an interesting problem to look at in that the combined effects of the salt 
and electric field induce cylindrical deformations [21] on such vesicles. 

The above discussions leads us, using equations (4) and (5), to write the Hamiltonian as 

K = \Jd 2 x{K t {V 2 hf+o t {Vh) 2 ). (9) 

This form is particularly useful for deriving the height-height correlation function, from which the 
roughness associated with the membrane can be deduced. 
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Figure 5: Variation of the lipid membrane's roughness versus the electric field for L = 40/im, where 
(To = lO -6 -^ in the (a) panel and do = 10 -7 -^ in the (b) panel. In each panel kq = 10 _20 (J),5 x 
10 _20 (J) and 10~ 19 (J) correspond to the top (blue), middle (red) and bottom (black) curves. 



3 The influence of an electric field on the roughness of lipid mem- 
branes 

The study of roughness in lipid membranes is facilitated by treating the evolution of such a system 
in terms of stochastic processes. One way to introduce stochastic behavior into the dynamics of the 
membrane is to treat the height function, h(x\,X2), as a stochastically fluctuating Wiener variable. 
This stochastic behavior can be communicated to the inclusions residing inside, and on the surface 
of the membrane. One may then proceed to obtain the height-height correlation function of the 
membrane in real space from which the roughness, associated with random changes in h, can be 
extracted. In the presence of an electric field and using the the Fourier representation of h, we 
arrived at the Hamiltonian [17] 

n= \lw? ^ Ktqi + atq2 } h{ci)h * {q) ' (10) 

where * denotes the complex conjugate. The static height-height correlation function can then be 
calculated as follows 

(MqjOMcrtO)) = ( 4 fc f 2 ) (2vr) 2 <5(q-q'). (11) 

V K tQ t otQ / 

The averaging is done with respect to the Boltzmann weight factor exp(— with k B being the 
Boltzmann constant. One can write the dynamic correlation function as [15] 

/i(q;t) = /i(q;0)e-^, 

leading to 

KB " (2vr) 2 J(q-q'), (12) 

where j(q) is the damping factor reflecting the long-range character of the hydrodynamic damping 
given by 

1 ' 3 



7(9) = ^ [^tq 6 + <J t q) . (13) 
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Figure 6: Variation of the lipid membrane's roughness versus the electric field for L = 500nm, 
where an = 10~ 6 ^ in the (a) panel and <7n = 10~ 7 ^j in the (b) panel. In each panel kq = 
1CT 20 (J),5 x 10~ 20 (J) and 10~ 19 (J) correspond to the top (blue), middle (red) and bottom (black) 
curves. 



Let us now suppose that the fluid surrounding the membrane has a viscosity n. It is then possible 
to write the Fourier transform of (12) in real space (i.e, (/i(x; t)h(x; 0))), representing the dynamical 
mean roughness of the membrane as 

o f f d 2 q d 2 q' p -x.(q-q') 



w = 1 1 (#(# (/l(q;t)r(q ' 0))e " x ' ' (14) 



Substituting from (12) and (13) we find 



2 1 fa , e ir > 



w =7T15 d Q 5 » ( 15 ) 

2TrpJj. K t qZ+o- t q 



where L is the linear size of the membrane and a is the molecular cut-off, of the order of a nanometer, 
and (3 = j^T- Here t represents the correlation time [15] and is of the order of 10 -4 s for the said 
values. Since there is no closed form expression for the integral given in equation (15), it should be 
computed by expanding the integrand as a series. The square root of w 2 is called root-mean-squared 
roughness (rms) when one assumes that (h(x,t)) = 0. Figure 5 and Fig. 6 demonstrate the variation 
of w with respect to an electric field. We have limited the electric field values to those presented 
in Table I. Here, the three different values for the bending rigidity kq = 10 _20 (J),5 x 10 _20 (J) and 
1CP 19 (J) correspond to the top, middle and bottom curves in each panel respectively. Note that in 
the (a) panels, Co = lO -6 ^- while in the (b) panels <to = lO -7 ^-. Also the typical values have been 
used as a = 2nm, T = 300A", and n = W' 3 ^. 

Note that the dimension of roughness in the figures is nanometer. The following observations 
worth mentioning: that the roughness values are of the order of a few hundred nanometers and that 
the decrease in kq and a® and the increase in L cause the roughness to increase. Here, it would 
be useful to mention that the observations of a typical membrane with L = 10 /im by AFM have 
yielded a roughness in the range 12-70 nanometers [6]. Also, results from molecular dynamics and 
monte carlo simulations of lipid membranes with sizes of the order of nanometers [12] show that they 
experience a roughness of the order of a few angstroms. For stronger fields, we see a dramatic change 
in the roughness, starting roughly at 30^ for an = lO -6 -^- and 10^ for an = 10~ 7 -t-. Therefore, the 
presence of an external electric field leads to the undulations or roughening of the membrane surface. 
Smaller values for correlation time lead to a lesser increase in rms roughness and vice versa. 
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4 Conclusions 

In this paper we have studied the behavior of lipid membranes in the presence of an external electric 
field. The negative contribution to the surface tension as a result of the application of the field, the 
effects of an external electric field on the roughening of membranes and an estimation of the thickness 
of charges aggregated on the membrane surface were reviewed. The dependence of bending rigidity 
on electric fields and its relation to the thickness mentioned above was studied. The increase in the 
membrane rms roughness for larger electric field values was calculated and shown to contribute to 
the roughening of the membrane surface. 
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